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We present the exact solution of Einstein's equation corresponding to a static and plane symmetric 
distribution of matter with constant positive density located below z — 0. This solution depends 
essentially on two constants: the density p and a parameter k. We show that this space-time finishes 
down below at an inner singularity at finite depth. We match this solution to the vacuum one and 
compute the external gravitational field in terms of slab's parameters. Depending on the value of «, 
these slabs can be attractive, repulsive or neutral. In the first case, the space-time also finishes up 
above at another singularity. In the other cases, they turn out to be semi-infinite and asymptotically 
flat when z — » oo. 

We also find solutions consisting of joining an attractive slab and a repulsive one, and two neutral 
O |' ones. We also discuss how to assemble a "gravitational capacitor" by inserting a slice of vacuum 

D ■ between two such slabs. 
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q-i I. INTRODUCTION 

U . 

bJQ, Due to the complexity of Einstein's field equations, one cannot find exact solutions except in spaces of rather high 
symmetry, but very often with no direct physical application. Nevertheless, exact solutions can give an idea of the 
qualitative features that could arise in General Relativity, and so, of possible properties of realistic solutions of the 
J> ■ field equations. 

\& | In this paper we want to illustrate some curious features of gravitation by means of a simple solution: the gravi- 
r**"* • tational field of a static plane symmetric relativistic perfect incompressible fluid with positive density located below 
z — 0. Because of the symmetry required, the exterior gravitational field turns out to be Taub's plane vacuum solution 
. The internal solution was also found by Taub Q. 

0^ ' Here we match both solutions.this corresponds to the plane symmetric counterpart of the Schwarzschild solution 
for a sphere of constant density 0. 

The solutions turn out to be attractive, repulsive or neutral depending on the value of a parameter. These space- 
times present some somehow astonishing properties without counterpart in Newtonian gravitation: 

Attractive solutions finish high above at an empty (free of matter) repelling boundary where space-time curvature 
diverges. These singularities are not the sources of the fields, but they arise owing to the attraction of distant matter 
as pointed out in [J|. 

Repulsive slabs explicitly show how negative but finite pressure can dominate the attraction of the matter. 
We also consider matching two different internal solutions. 

In Sec. II we present a simple and complete derivation of Taub's interior solution. In Sec. Ill we make a detailed 
study of it. In Sec. IV we discuss how solutions can be matched. 

Throughout this paper, we adopt the convention in which the space-time metric has signature (— + + +), the 
system of units in which the speed of light c = 1, Newton's gravitational constant G — 1 and g denotes gravitational 
field and not the determinant of the metric. 



II. THE TAUB INTERIOR SOLUTION 



In this section we consider the solution of Einstein's equation corresponding to a static and plane symmetric 
distribution of matter with constant positive density and plane symmetry. That is, it must be invariant under 
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translations in the plane and under rotations around its normal. The matter we shall consider is a perfect fluid 
satisfying the equation of state p = p , where po is an arbitrary positive constant. The stress-energy tensor is 

T a b = (P + P) U a Ub + P g a b , (1) 

where u a is the velocity of fluid elements. 

Due to the plane symmetry and staticity, following [l| we can find coordinates (t, x, y, z) such that 

ds 2 = -G(z) 2 dt 2 + e 2V ^ (dx 2 + dy 2 ) + dz 2 . (2) 

That is, it is the more general metric admitting the Killing vectors d x , d y , xd y — yd x and dt- 
The non identically vanishing components of the Einstein tensor are 

G tt = -G 2 (2 V" + 3V' 2 ) (3) 

Gxx = G yy = e 2V (g"/g + g'/g v + v" + v 2 ) , (4) 

g zz = V (2 g'/g + V) , (5) 

where a prime (') denotes differentiation with respect to z. 

On the other hand, since the fluid must be static, u a = (g, 0, 0, 0), so 

T ab = di ag (pg 2 ,pe 2V ,pe 2V ,p) , (6) 

where p depends only on the z-coordinante. Thus, Einstein's equations, i.e., G a b — 8nT a b, are 

2V" + 3V' 2 = -87rp, (7) 

g"/g + g'/g v + v + v 2 = 8n P , (8) 

v'(2g'/g + v') = 87rp. (9) 

Moreover, V ' a T ab = yields 

P ' = -( P +p)g'/g. (io) 

Of course, due to Bianchi's identities equations, ([7]), ([5]), © and (JTUJ) are not independent, so we shall here use only 
©, ©, and (dHD- 

Since p is constant, from (|10[) we readily find 

P = c p /g(z)~p, (ii) 

where C p is an arbitrary constant. 

By setting W(z) — e 3V ^ z ^ 2 , we can write (JTJ) as W" = —6npW, and its general solution can be written as 

W(z) =Ci sini-s/fhrp z + C 2 ), (12) 
where C\ and C2 are arbitrary constants. Therefore, we have 

V(z) = In (C x sin {y/f^p z + C 2 )) 3 . (13) 
Now, by replacing (fTTj) into |9]) we get the first order linear differential equation which g(z) obeys 



= —\J 67TO [tan u H — cot u ) g + J 6np — tan it, (15) 
V 3 j p 

where u — ^6np Z + C2, and in the last step we have made use of (fl~3|) . The general solution of (fT"4l) can be written as 

cos it / C p f (sinu) 4 / 3 

y — T- \TR const. H / — ^-du 

smii ''^ \ p J (cosu)^ 



G 



cos it 3Cr, .9 „ / 2 13 



(sinit) 1/3 7p 



f- sin 2 u 2 F 1 (l,-;— ; sin 2 u), (16) 
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where C3 is another arbitrary constant, and 2-^1(0, 6; c; z) is the Gauss hypergeometric function (see for example 0]). 
Therefore, the line element ^ becomes 

ds 2 = -G{zf dt 2 + (Ci sinu)^ (dx 2 + dy 2 ) + dz 2 , (17) 

where Q(z) is given in (|16[) and u = i/Enp z + G%. Thus, the solution contains five arbitrary constants: p, C p , C\, C2, 
and C3. The range of the coordinate z depends on the value of these constants. 

For p > 0, C p > and C3 > 0, this solution was found by Taub Q. When C p — 0, it is clear from (jTT[) that 
p(z) = —p, and the solution (|17|) turns out to be a vacuum solution with a cosmological constant A = 8irp [5|. For 
C p = and p — > 0, by an appropriate choice of the constants we can readily see that (|17p becomes 

ds 2 = -(1 - Zgz)~i dt 2 + (1 - 3gz)3 (dx 2 + dy 2 ) + dz 2 , 

—00 < t < 00, — 00 < x < 00, —00 < y < 00, < 1 — 3gz < 00 , (18) 

where g is an arbitrary constant. In (|18p the coordinates have been chosen in such a way that it describes a homoge- 
neous gravitational fieldo pointing in the negative ^-direction in a neighborhood of z = 0. The metric (|18p is Taubs's 
vacuum plane solution [l[ in the coordinates used in 0]. 

Nevertheless, the solution (|17[) has a wider range of validity. For abnormal matter some interesting solutions also 
arise, but the complete analysis turns out to be somehow involved. So, for the sake of clarity, we shall confine our 
attention to positive values of p and C p ^ 0, leaving the complete study to a forthcoming publication Q. 

Notice that the metric (|17[) has a space-time curvature singularity when u = 0, since straightforward computation 
of the scalar quadratic in the Riemann tensor yields 

RabcdR ahcd = 4 (g" 2 + 2 g' 2 v 12 ) /g 2 + 4 (2 v" 2 + 4 v"v' 2 + 3 v' 4 ) 

+ + + (!_,)), (M , 

so R a bcdR abcd — ► 00 when m — > 0. 



III. THE PROPERTIES OF THE FUNCTION Q{£) 

In this section we shall study in detail the properties of the solution for the case p > and C p ^ 0. 

Now, it is clear from ([7]), ([5]), ([5]) and (TIT))) that field equations are invariant under the transformation z — > 
±z + zo, i.e., z-translations and mirror reflections across any plane z=const. Thus, if {g(z) 1 V{z),p{z)} is a solution 
{g(±z + zo), y(±z + zo),p(±z + zo)} is another one, where zo is an arbitrary constant. Therefore, taking into account 
that u — y/6irp z + C2, without loss of generality the consideration of the case < u < ir/2 shall suffice. 

By an appropriate rescaling of the coordinates {t,x,y}, without loss of generality, we can write the metric (117[) as 

ds 2 = -g(z) 2 dt 2 + sins u (dx 2 + dy 2 ) + dz 2 , 
— 00 < t < 00, —00 < x < 00 , — 00 < y < 00, < u = \J &np z + C2 < 7r/2, (20) 

and (fl!)j) as 

, kC v cosu 3Cp . 9 „ / 2 13 . 9 \ ,„ . 

G 2 = — — TJ ^ + ^ L sin 2 M 2 iq(l,- ; — ; sin 2 u), 21 
p sin 1/J !i 7p V 3 6 / 

where k is an arbitrary constant. 

By replacing (f2"Tj) into (TTTJ) , we see that the pressure is independent of C p . On the other hand, since G(z) appears 
squared in the metric, it suffices to consider C p > 0. Furthermore, rescaling again the coordinate t, we may set 
C p = p. Thus, (j2T|) becomes 

_ , . _ , . cosu 3.9 „/ 2 13 . , \ 

Q(z) = G k (u) = «l : +-sm 2 « 2 Fi 1,-;— ;sin 2 tx , 22 

sin 1/d tt 7 V 3 6 J 

where G K (u) was defined for future use, and we recall that u = y/WFpz + C2. Furthermore, pT|) becomes 



p(2) = p (1/5(2) -1). 



(23) 







tt/4 

(a) K > &crit 



tt/4 

(b) ft = t^crit 




Therefore the solution depends on two essential parameters, p and n. We shall discuss in detail the properties of 
the functions Q{z) and p(z) depending on the value of the constant k. 

Now, the hypergeometric function in the last equation is a monotonically increasing continuous positive function 
of u for < u < 7r/2, since c — a — 6 = 1/2 >0. Furthermore [|| 

^1,§ ; ^0)-1, and^l,?;^),!. (24) 

Therefore, we readily see from (|2"21 that, no matter what the value of n is, g(z)\ u —„/2 — 1 an d we get then from (|2"5|) 
that vanishes at tt = 7r/2. 
On the other hand, since 

G(u) = ku- 1 / 3 + 0(u 5 / 3 ) as u^O, (25) 

£(-z)|u=o = if k = 0, whereas it diverges if k ^ 0. 

For the sake of clarity, we shall analyze separately the cases n > 0, k — 0, and K < 0. 



i. « > o 



In this case, it is clear from (|22|) that is positive definite when < u < ir/2. On the other hand, from {8]) and 
)]) we get 

g" = g'v'- gv" = - (v" + + ^p) Q + ^c p 

= V' 2 g + 4irC p , (26) 
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where we have made use of (| 14[) and |7]). Then also Q" is positive definite inO < u < n/2, and so Q' is a monotonically 
increasing continuous function of u in this interval. 

Now, taking into account that Q' = d z Q — ^6np d u Q, a straightforward computation from (f22f shows that 



= _^V^ u -4/3 + 0(u 2/3j asw ^Q, (27) 

and 

^'(z)|n=7r/2 = \f^P{^crit ~ «) , (28) 

where 

Kcrit = %Ar(7/6)/r(2/3) = 1.2143 . . . (29) 

If K > K cr i tl G' is negative for small enough values of u and non-positive at u = 7r/2, hence £7' is negative in 
< u < 7r/2, so C/(z) is decreasing, and then Q{z) > G(z)\ u=7r /2 = 1 in this interval (see FigQTa) and FigQTb)). 

For K cr u > k > 0, Q' is negative for sufficiently small values of u and positive at 7r/2. So there is one (and only one) 
value u m where it vanishes. Clearly Q{z) attains a local minimum there. Hence, there is one (and only one) value uq 
(0 < it < 7r/2) such that Q{z)\ u=Uo = G(z)\ u=7T /2 — 1, and then Q(z) < 1 when u < u < n/2 (see Figfjjc)). 

Since Q{z) > 0, it is clear from (|23|) that p(z) > if Q(z) < 1, and p(z) reaches a maximum when Q{z) attains a 
minimum. 

Therefore, for k > n cr i t , p(z) is negative when < u < n/2 and increases monotonically from — p to (see FigOJa) 
and FigETTb)). 

On the other hand, for n crit > k > 0, p(z) grows from — p to a maximum positive value when u — u m where it 
starts to decrease and vanishes at u = n/2. Thus, p{z) is negative when < u < u and positive when uq < u < n/2 
(see FiglUc)). 

2. k = 

In this case, it is clear from (|22|) that Q monotonically increases with u from to G(z)\ u =tt/2 = 1- Therefore, p is 
a monotonically decreasing positive continuous function of u in < u < n/2 (see Fig[TJd)). Furthermore, at u = it 
diverges, since 

p(z) - ^-u~ 2 -> +c» as u — >• 0. (30) 
5. K < 

In this case, we see from (|27p that CJ' is positive when u takes small enough values, and from l|28p that it is also 
positive when u is near to n/2. 

Now, suppose that Q'{z) attains a local minimum when u = u\ (0 < u\ < n/2), then G" (z)\ u=Ul = 0. Hence, we 
get from (|2"rJ)) that G(z)\ u=Ul < 0, and taking into account that V'{z)\ u=Ul — 2y/6np/3 cot u\ > 0, we see from (TH]) 
that ^'(^)|„_ Ul > 0. Thus, we have shown that G'{z) is a continuous positive definite function when < u < n/2 if 
k < 0. 

Therefore, in this case, G{z) is a continuous function monotonically increasing with u when < u < n/2. Since it 
is negative for sufficiently small values of u and 1 when u = n/2 it must vanish at a unique value of z when u = u K 
(say). Furthermore G{z) < 1 when < u < n/2. 

Clearly, we get from (|22p that u K is given implicitly in terms of ac through 

3 sin 7 / 3 u K / 2 13 a \ 

« = -= 2^1 1,-;— ;sm m k . (31) 

7 cosu K V 3 6 / 

We can readily see from (f3"Tj) that u K is a monotonically decreasing function offtin— oo<k<0, and it tends to n/2 
when At — - > — oo and to when k — > _ . 

From (|23|). it is clear that p(z) diverges when u — u K . Furthermore, |23|) also shows that p(z) < when G(z) < 0, 
and taking into account that G(z) < 1, that p{z) > when G(z) > 0. Therefore, p{z) is negative when < u < u K 
and positive when u K < u < n/2 (see FigJ^J). 



6 




On the other hand, we see from (|19p that, when k is negative, another space-time curvature singularity arises at 
u K (besides the one at u = 0) since p diverges there. 

Therefore, if k is negative, the metric (I2U1) describes two very different space-times: 

(a) For < u < u K , the whole space-time is trapped between two singularities separated by a finite distance 
^]%-Kpu K . This is a space-time full of a fluid with constant positive density p and negative pressure p monotonically 
decreasing with u, and p(z)\ u —o = —p and p(z) — > — oo as u — > u K . 

(b) For u K < u < 7r/2, the pressure is positive and monotonically decreasing with u, p{z) -too as a -> ii K and 
p(z)\u=7r/2 = 0. 



IV. MATCHING SOLUTIONS 



In this section, we shall discuss matching the interior solution to the vacuum one, as well as joining two interior 
solutions facing each other. 

Since the equations are invariant under z-translation, we can choose to match the solutions at z = without losing 
generality. So we select C2 = 7r/2, and then (|22|) becomes 

Q{z) = G K (^Tpz + tt/2) = —k S1 ^^f } r + I cos 2 (V6^z) 2 Fi(l ,^ H ; cos 2 (v^ Z )) . (32) 

cos 1 ' d (yo7rp z) 7 V o 6 / 

Therefore, the metric (|2T)|) reads 

ds 2 = -G K (^67rpz + tt/2) 2 dt 2 + cos^ (^6tt p z) (dx 2 + dy 2 ) +dz 2 , 
—00 < t < 00, —00 < x < 00, -co < y < 00, — V7r/24p < z < . (33) 

Notice that it holds that g u (0) = -£(0) 2 = -1, g xx (0) = g yy (0) = 1, and from {23]) we get p(0) = 0. 
Furthermore, we shall impose the continuity of d z gu, and from (f2"5| we have 

d z gtMVnt en or = -2^(0)^(0) = -2 v /6^(k c „ 4 - «) . (34) 



A. Matching solutions to vacuum solutions 



Here, we shall consider the gravitational field of a planar symmetric distribution of matter with constant density p 
sitting below z = 0. This is the plane symmetric counterpart of Schwarzschild's solution for a sphere of incompressible 
fluid [|. 
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In this case, we see from (IT8l) that the corresponding exterior solution for z > is 

ds 2 = -(1 - 3gz)- 2 ^ 3 dt 2 + (1 - 3. 9 z) 4 / 3 (dx 2 + dy 2 ) + dz 1 



-oo < i < oo, — oo < a; < oo, — oo < y < oo, < z < 



1/35 if 5 > 
oo if g < ' 



(35) 



which for g > is the region finishing at the singularity of Taubs's vacuum plane solution, whereas it is the asymptotic 
flat tail of it for g < 0. It describes a homogeneous gravitacional field — g in the vertical (i.e., z) direction [4j. 

Since gtt(0)\exterior = -1, gxx(0)\exterior = g yv (0)\exterio7- = 1 the continuity of the metric is assured, and concerning 
the derivative we have 



9 z gtt (z) | exterior = ~2g (1 " 3.gz) 5/3 . 

Therefore, by comparing with (|34| . we see that the continuity of d z gu yields 



g = y/6np (k c 



k) , 



(36) 



(37) 



which relates the external gravitational field g with matter density p. 

Thus, we can readily see from (|37|) that the slab is attractive if k < n C rit, and the whole space-time is trapped 
between two singularities: the inner one discussed in the previous section and the outer one at z = l/3<?. 

If k < 0, the range of the z-coordinate is — (tt/2 — u K ) / y/Qup < z < l/3g, where u K (0 < u K < ir/2) is given 
implicitly in terms of k through (l3T|) . Inside the slab, the pressure is always positive, and it diverges deep below at 
the inner singularity (see Fig(2J) . Its depth is 



(tt/2 - w K )/ v /67rp. 



By using (|31|) . we can write k in terms of d 

3 cos 7 / ' 3 (y/6^d) 
7 sm(^/6np d) 



('■ 



2 13 

3 ; y ; 



\^femp~d) 



(38) 



(39) 



From (|3"T|) . (|2^|) and (|3^1) . we can write the external gravitational field g in terms of the matter density p and the 
depth of the slab d 



,' V^r(7/6) 3 cos 7 / 3 U/6^d) _/ 2 13 2/ /«— ^ 



(40) 



If < k < K cr it, the z-range is —\J^j2Ap < z < l/3g and inside the slab, the pressure is always finite, but it is 
negative deep below and p = —p at the inner singularity (see Fig[T{c)). 

Therefore, when k < n C rit, we see how the attraction of the distant matter lying below z — shrinks the space-time 
in such a way that it finishes high above at the outer empty singular boundary (at z — l/3.g). 

It can readily be seen from (|3T[) that, if n > K C rit, g is negative and the slab turns out to be repulsive. In this case, 
the space-time is semi-infinite and asymptotically flat when z — > oo (see j4j). Inside the slab, the pressure is always 
finite and negative, and p = —p at the inner singularity (see FigfTJa)). 

If « = K cr it it is gravitationally neutral, and the exterior is one half of Minkowski's space-time. 

Two remarks are in order. First, notice that the maximum depth that a slab with constant density p can reach 
is y/ir/2Ap : being the counterpart of the well-known bound M < 4R/9 (R < l/^/Sirp) which holds for spherical 
symmetry. 

On the other hand, notice that the derivative of g X x(z) = g y y(z) has a discontinuity at the surface (z = 0) as it 
occurs with d r g rr {r) in the Schwarzschild's case 0, since 



d z g X x(0)\i 



®z9yy (0)1 inte 



4-y / Tj7rp 



cos 3 (y/ 6irpz) sin(\/ 6irp z)|z=o = . 



(41) 



but 



(0)1 exterior — $zQyy (0)1 exterior — 

-4g(l — 3gz) ' | z=0 



(42) 
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B. Matching two slabs 

Now we consider two incompressible fluids joined at z = where the pressure vanishes, the lower one having a 
density p and the upper having a density p' . Thus, the lower solution is given by (|33[) . By means of the transformation 
z — > —z, p — > p' and n — > «/ we get the upper one 

ds 2 = -G K ,(n/2- ^6Trp'z) 2 dt 2 + cos^ ( ^6^' z) (dx 2 + dy 2 ) + dz 2 , 
— oo <t< oo, — oo < x < oo, — oo < y < oo, < z < y 7r/24// . (43) 

From (|33[) and (|43[) . we can readily see that g u {z), g X x{z) and d z g xx {z) are continuous at z = 0. Furthermore, from 
([34]) we see that the continuity of d z gtt requires 

y/p(Kcrit ~ «) = -VV («crit ~ «') • (44) 

Thus, if one solution has a k greater than K cr n , the other must have it smaller than « cr jt. Therefore, the joining is 
only possible between an attractive solution and a repulsive one, or between two neutral ones. 

It is easy to see that we can also insert a slice of arbitrary thickness of the vacuum solution (| 1 8|) between them, 
obtaining a full relativistic plane "gravitational capacitor" . For example, we can trap a slice of Minkowski's space-time 
between two solutions with k = n cr i t . 



V. CONCLUDING REMARKS 



We have done a detailed study of the exact solution of Einstein's equations corresponding to a static and plane 
symmetric distribution of matter with constant positive density. By matching this internal solution to the vacuum 
one, we showed that different situations arise depending on the value of a parameter k. These simple space-times 
turn out to present some somehow astonishing properties. 

For k < K cr it, the attraction of the distant matter shrinks the space-time in such a way that it finishes high above 
at an empty singular boundary as pointed out in [J] . This space-time also finishes down below at another singularity. 

For k < 0, we have explicitly computed the external gravitational field g in terms of the density p and the sickness 
d of the slab. 

For k > K cr it, negative but finite pressure dominates the attraction of the matter and the slab turns out to be 
repulsive. 

We showed that the maximal sickness that these slabs can have is y / 7r/24 / o. 

We have also discussed matching an attractive slab to a repulsive one, and two neutral ones. We also comment 
on how to assemble relativistic gravitational capacitors consisting of a slice of vacuum trapped between two of such 
slabs. 
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